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ABSTRACT 

We study the imprints on the formation and evolution of cosmic structures of a particular class 
of dynamical dark energy models, characterized by an oscillating equation of state. This in- 
vestigation complements earlier work on the topic, that focused exclusively on the expansion 
history of the Universe for such models. Oscillating dark energy cosmologies were introduced 
in an attempt to solve the coincidence problem, since in the course of cosmic history matter 
and dark energy would have had periodically comparable energy densities. In this class of 
models the redshift evolution of the equation of state parameter w{z) for dark energy is char- 
acterized by two parameters, describing the amplitude and the frequency of the oscillations 
(the phase is usually set by the boundary condition that w(z) should be close to —1 at re- 
cent times). We consider six different oscillating dark energy models, each characterized by 
a different set of parameter values. For one of these models w{z) is lower than —1 at present 
and larger than — 1 in the past, in agreement with some marginal evidence from recent type 
la supernova studies. Under the common assumption that dark energy is not clustering on the 
scales of interest, we study different aspects of cosmic structure formation. In particular, we 
self-consistently solve the spherical collapse problem based on the Newtonian hydrodynam- 
ical approach, and compute the resulting spherical overdensity as a function of cosmic time. 
We then estimate the behavior of several cosmological observables, such as the linear growth 
factor, the Integrated Sachs-Wolfe (ISW) effect, the number counts of massive structures, and 
the matter and cosmic shear power spectra. We show that, independently of the amplitude and 
the frequency of the dark energy oscillations, none of the aforementioned observables show 
an oscillating behavior as a function of redshift. This is a consequence of the said observables' 
being integrals over some functions of the expansion rate over cosmic history, thus smooth- 
ing any oscillatory features in w{z) below detectability. We also notice that deviations with 
respect to the expectations for a fiducial ACDM cosmology are generically small, and in the 
majority of the cases distinguishing an oscillating dark energy model would be difficult. Ex- 
ceptions to this conclusion are provided by the cosmic shear power spectrum, which for some 
of the models shows a difference at the level of ^ 10% over a wide range of angular scales, 
and the abundance of galaxy clusters, which is modified at the ~ 10 — 20% level at z > 0.6 
for future wide weak lensing surveys. 
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1 INTRODUCTION 

In recent years an increasingly large body of observations con- 
firmed the general framework of a standard cosmological model 
based on General Relativity. Accordingly, right after the Big Bang 
the Unive rse experienced an accelerated expansion phase, dubbed 
inflation ( iGutrill98ll : [Lindelll982l: IZelnikovllT99ll) . during which 
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quantum fluctuations were amplified to produce tiny perturbations 
in the matter distribution, whose imprints can nowadays be ob- 
served in the Cosmic Microwave Background (CMB) temperature 
map. Later on, and due to gravitational instability, these seed fluc- 
tuations grew up, giving rise to the web of cosmic structures that we 
observe to day. After inflation, the Universe e xperienced a period of 
reheating jShtanov. Traschen & Brandenber ger 1995; Kaiser 199^ 
with the formation of light elements. Current observations of the 
CMB and of the luminosity distances of type la Supernovae (SNe 



2 F. Pace et al. 



la) show that the geometry of the Universe is spatially flat, in ac- 
cordance with the predictions of the inflationary paradigm, and fur- 
thermore showing that the Universe is currently undergoing another 
accelerated expansion phase. 

After the first detection of an accelerating expansion rate 
at low redshift (p oss ibly z < 0.5, IShap iro & T urned l2006h by 



iRiess et alj h998l) and lPerlmutter et al.l (^99 ) many other different 
and independent studies led to the same conclusions making this in- 
ference very solid. In partic ular, evidence for an acce lerated expan- 
sion comes from the CMB jjaffe et al.ll200ll : lKom"ats u et al. 2011) 
and the Integrated Sachs-Wolfe (ISW) effect ( IHo et alJi2008ll . the 
Large Scale S tructure (LSS) and the Baryon Acousti c Oscilla- 
tion (BAO) (Eisensteinet al. 2005: Perci val et alj |2010). globular 
clusters (iKraus s &. Chabover 2003), galaxy clus ters CHaiman et al.l 
l200ll;lAllenet al. 2004, 2008; Wan g et alj|2004) and weak lensing 
jHoekstra et"ar.ll2006i : Jarvis et al.ii200^ . 

As a homogeneous and isotropic General Relativistic model 
universe filled with matter is unable to reproduce the observed ac- 
celerated expansion, three different explanations have been pro- 
posed to account for it. One possibility consists in putting aside the 
hypothesis of homogeneity on large scales: these models are de- 
scribed by the Lemaitre-Tolm an-Bondi (UTB) m etric or are based 
on the idea of backreaction jKolbetal.ll200^ . A second possi- 
bility is to suppose that on very large scales General Relativity 
breaks down and gravity is modified. In this case we will be in the 
need for a new theory of gravity and General Relativity would be 
only the small-scale limit of a more profo und theory. Examples of 
this id ea are the f{R) models ( Amendol a et alj|2007l:lstarobinskvl 



l2007h . brane models ('Oeffayet 2001; Dvali e t al.ll20q(]h. and the 



f{T) models {B engochea & Ferraro 2009; Lin derl2010l : lDent et al.1 
12011 ': Zheng & Huang|l2011l) . Finally, one could assume that Gen- 
eral Relativity is correct but the low-z Universe is dominated by 
some kind of exotic fluid with negative pressure, the dark energy. 
Specifically, if dark energy constitutes ~ 70% of the matter-energy 
content of the Universe, from the second Friedmann equation it 
turns out that its equation of state parameter w would need to be 
w < —1/2 in order to provide accelerated expansion. 

In the concordance cosmological model the role of dark en- 
ergy is played by the cosmological costant A, having a redshift- 
independent equation of state parameter w — —1 and commonly 
interpreted as the energy density of the vacuum. Even though this 
A-Cold Dark Matter (ACDM henceforth) model is now the stan- 
dard reference framework in cosmology, it suffers from some fun- 
damental theoretical problems, that can be summarized by the fol- 
lowing questions. 

• Why is the energy density implied by the cosmological con- 
stant much smaller than the theoretically expected vacuum energy 
density? 

• Why is the dark energy density comparable to the dark matter 
density only today? 

The last one is also known as the coincidence problem. In order 
to solve or at least alleviate these issues, it is possible to iden- 
tify dark energy with the energy density of a minimally coupled 
scalar field (named quintessence), that evolves through cosmic 
time as dictated by its own potential. This gives rise to a redshift- 
dependent equation of state parameter w{z), hence making at least 
the coincidence problem less severe. These dynamical dark energy 
mo dels can be roughly gr ouped into two classe s, tracking mod- 
els (ISteinhaj-dtet alJll999l) and sc aling models ( iHalli well! [19871 : 
IWands et alJl993l : IWettericlJl995l) . 

Models with an oscillating equation of state were introduced 



to solve the coincidence problem, because the present accelerated 
expansion phase would just be one of the many such phases oc- 
curring over cosmic history, especially at early times. Moreover, 
oscillations would more naturally accommodate the crossing of the 
phantom barr ier, w = —1 as it is marginally suggested by recent 
observations ^A lam et tdj|2004l: lAUen et al.ll2004 loic us & Repkd 
[2004; Riess et a lJl2004lFen^ et al Jl2005l : lHuterer & C ooray 200i 
■ Choudhurv & Pa dmanabhan 2005). In the framework of particle 
physics, it is possible to have an oscillating quintessence poten- 
tial if one considers a pseudo-Nambu Goldstone boson field when 
it has rolled through the minimum. As models for dark energy, os- 
cillating scalar fields were propo sed byjDutta & Scherrer (200^ 
Ijohnson & Kamionkowski' (2008), and Gu| j2008l) . iLazkoz et 
1 20051) and lKurek et al. (2008) found a better agreement with SNe 
la data if an oscillating equation of state is used instead of the cos- 
mological constant or an equation of state linearly dependent on 
the scale factor. Further indi cations in the sam e direction come 
from the study performed by iRiess et ^ l l2007h . Previous works 
on this topic focused mainly on the expansion history of the Uni- 
verse and marginally on the linear perturbation theory in the frame 
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work of oscillating quintessence (Feng et al. 2005; Xia et al. 200: 
Barenboim & Lykken 2006; Barenboim et al. 2006; Kureket 
2008il201(]| : lLan et all20IQl . 



The novelty of this work lies in the fact that we explore sig- 
natures of an oscillating equation of state w{z) in the (non-linear) 
growth of cosmic structures, thus extending and complementing the 
majority of foregoing studies. The main idea we explore here is to 
find out cosmological observables based on structure formation that 
can hint toward oscillating quintessence even though the expansion 
of the homogeneous and isotropic background does not. The rest of 
this paper is hence organized as follows. In Section |2] we describe 
the formalism of the spherical collapse used to derive important pa- 
rameters for the formation and evolution of structures. In Section[3] 
we describe and motivate several parametrizations used in order to 
describe the oscillating dark energy and in Section |4] we present 
results for the different observables we considered. Section[5]is de- 
voted to our conclusions. In Appendix lAl we present implementa- 
tion details of the code used to evaluate the linear growth factor and 
the evolution of the spherical overdensity. 



2 SPHERICAL COLLAPSE MODEL 

Despite its simplifying nature, the model describing the collapse 
of a uniform non-rotating spherical overdensity in a cosmological 
setting provides numerous insights on the actual process of struc- 
ture formation. For instance, the linear density contrast extrapo- 
lated at the spherical collapse time provides a fair approximation 
for the threshold at which actual perturbations can collapse to form 
bound structures. Thus, in this Section we sketch the derivation of 
the relevant equations for the spherical collapse model under the 
assumption that only dark matter can form clumps, while dark en- 
ergy is just present as background fluid. For fur ther details w e re- 
fer to the curren t literature on the topic (see, e.g.iBemardeaj|l994l: 



Ohta et al.l2003l,l20()4l:lMota & van de Bruckl2004lNunes & Mot3 



200d: lAbramo et alj2007l : |Pace et al.ll2010l) . 

Rather than studying the time evolution of the radius of the 
collapsing sphere, we study directly the time evolution of the over- 
density. This procedure proves to be numerically more stable and 
less prone to errors than the classical approach based on the radius 
evolution jPace et alllloTol) . We consider a perfect fluid described 
by the energy-momentum tensor T''", satisfying the local conser- 
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vation laws expressed by V^T'"^ = 0. This set of four equations 
encapsulates both the continuity and the Euler equations, while 
from Einstein's field equations it is possible to derive a relativis- 
tic generalization of the Poisson equation. In a more explicit form 
these expressions read 

■{pv) + ^V ■v = Q, (1) 



the non-linear evolution of the density perturbation up to the time 
of virialization, we need to evaluate the tum-around scale factor 
ttta, defined as the time when the radius of the sphere reaches its 
maximum, detaches from the overall expansion of the Universe, 
and collapses afterwards. Using then the virial theorem and energy 
conservation considerations, the vi rial overdensity Ay ca n be de- 
rived according to the discussion of lMaor & L^iavl j2005t) . 



^ + {v- v)i;+ v$ = . 
at 
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. 



(2) 



(3) 



where p, P, v and <1? are the density, the pressure, the velocity and 
the gravitational potential of the fluid. 

For the average background matter density the following continuity 
equation holds, 



p + 3H{p+-] =0, 



(4) 



where p = 3//^f2m/87rG is the background matter density, H 
is the Hubble parameter, and Qm is the matter density parameter. 
Since for ordinary matter and dark matter the pressure contribution 
is negligible, from now on we will set P = 0. 

Assuming spherical symmetry and perturbing the physical 
quantities appearing in the previous set of equations (density, veloc- 
ity and gravitational potential) around their background values, we 
obtain the following exact non-linear differential equation, describ- 
ing the evolution of the matter density perturbation (5 as a function 
of the cosmic time. 



. 4 ^2 



-4ttGp5{1 + S) = 0, 



(5) 



3 1 + 5 

We stress that Eq. l[5) is valid also for large density contrasts, deep 
in the non-linear regime, as long as spherical symmetry is satisfied. 
By restricting to 5 <C 1 instead, at first order Eq. ^ reads 



S + 2HS - 47rGp S^O, 



(6) 



and it coincides with the differential equation commonly used to 
determine the linear growth factor^ 

As explained in detail in IPace et al in order to de- 

termine the linear density perturbation threshold for spherical col- 
lapse, (5c, one should solve Eq. ^ with suitable initial conditions, 
namely the initial overdensity and velocity of the perturbation. In 
order to find the initial overdensity we take into account that at the 
time of the collapse of the object all the matter is concentrated in 
one point, therefore formally S — )■ +oo. Hence by fixing the time, 
or scale factor, of collapse ac, with a root-search method it is pos- 
sible to determine the initial overdensity Si such that the solution of 
the non-linear Eq. (O diverges at ac . Once the initial overdensity is 
found, we compute the initial velocity as detailed in Appendix IaI 
and use them both as initial conditions for the linear Eq. (|6}. When 
integrated up to Oc the latter returns us the linear density contrast 
corresponding to the time of spherical collapse, 5c . 

In order to determine the virial overdensity Av, representing 



3 OUTLINE OF THE COSMOLOGICAL MODELS 

As outlined in Section[T] at the moment there is no explanation for 
dark energy in terms of fundamental physics, therefore all the mod- 
els that we explored in the present work are purely phenomenolog- 
ical and the values of their parameters are generically adjusted such 
that certain classes of cosmological observables (most commonly 
the luminosity distance of SNe la and the CMB temperature power 
spectrum) are well reproduced. These models are described in the 
present Section. As a fiducial reference cosmology we assume the 
standard flat ACDM model. The cosmological parameters are set 
to fim.o = 0.27 4, Q.Q,o = 0.726 and h = 0.7 in accorda nce with 
WMAP-7 data jKomatsu et al.ll201 ll : Earson et alj|201 ll) and the 
Supernova Legacy Survey 3 data (see Sullivan et al. 20111) . These 
same parameter values are also kept intact for all the dynamical 
dark energy cosmologies investigated in this paper. 

The amplitude of the primordial matter power spectrum in the 
fiducial ACDM cosmology is selected in order to attain a given 
value of the quadratic deviation on a comoving scale of 8 Mpc//i, 
(Tg = 0.8. In all the other dynamical dark energy models that we 
considered the normalization is scaled according to 



0"8,DE = 



,{z = 0) 



0"8,ACDM , 



(7) 



(Sc,acdm(.z — 0) 

where 5c is the linear overdensity parameter extrapolated at spher- 
ical collapse (see Section |2] above). In this way, the exponential 
tail of the dark matter halo mass function at redshift zero is con- 
served, and hence the abundance of massive structures at present 
times, which is arguably well defined from the observational point 
of view, is the same for all models. We show the values for the 
normalization of the different models in Table [T] 

In the present work we analyzed six different dark energy cos- 
mologies with an oscillating equation of state parameter w{z). In 
the first five of them w{z) has the same functional form but dif- 
ferent values of the free parameters. The sixth model has instead a 
different functional form for w{z), although still presenting oscil- 
lations. The functional form for the first fiv e cosmologies is tLinde3 
l200d lLazkozetalj2010l : iFeng et al.ll2006h 



w{a) = Wo — A sin (B In a + 9) 



(8) 



where a = 1/(1 + 2) is the scale factor, A determines the ampli- 
tude of the oscillations, B gives their frequency while S is a phase 
shift. As can be easily seen, the value of the equation of state pa- 
rameter today is 'w{a = 1) — wo — j4sin(6'), which equals wq 
if the phase is ^ = 0. Model six is meant to generalize the CPL 
parametrization jChevallier & Polarski|[200ll : iTinder 2003) m or- 
der to avoid the future unphysical divergence of the dark energy 
equation of state typical of this model. In this case, the function 
w{z) can be written as 



w{a) =wo- A{aB sin (l/a) + 6*) 



(9) 
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Table 1. Values of the free parameters for the dark energy equation of state 
and for the matter power spectrum normalization erg. The numbers in th e 
last c olumn conspond to the following references: (1) : Linder 120 061). 
(2) : iLazkoz et alj i2010l) . (3) : iFeng et alj l2006l) . (4) : .Ma & ZhmS 
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-0.9 


0.07 


5.72 


0.0 


0.7989 


(1) 


2 


-0.9 


0.07 


2.86 


0.0 


0.7986 


(1) 


3 


-0.9 


0.15 


1.0 


0.0 


0.7983 


(1) 


4 


0.0 


1.0 


0.06 


7r/2 


0.7999 


(2) 


5 


-1.0 


1.5 


0.032 


57r/18 


0.8012 


(3) 


6 


-1.061 


0.041 


1.0 


— sin(l) 


0.8001 


(4) 



Going into the distant future, we have a — > oo and taking the limit 
of the previous equation we observe that it asymptotes to 

U)(a) = mo - A{B + 9) , (10) 

thus rem aining finite. For mo re details on this model we refer the 
reader to lMa & Zhan3 j201lh . 

In Table [T] we summarize the values for the free parameters 
characterizing each cosmology. We also quote the paper where the 
expansion history of the Universe resulting from that specific dark 
energy model has been studied in detail. The quoted values are the 
best fit to certain classes of cosmological observables considered 
by those authors. In the literature other parametrizations of w{z) 
showing an oscillatory behavior can be found (see for example 
iKurek et al. 2008 . 2010). However they represent only local fits to 
the expansion history, giving unphysical divergences in the distant 
past. For this reason we did not include them in our analysis. 

In Figure [T] we show the redshift evolution of the equation of 
state parameter 'w{z) for the different models studied in this work. 
We refer to the caption for the description of the line styles adopted. 
As it appears evident from the Figure, models 1 and 2 have the same 
amplitude but the frequency of oscillations changes by a factor of 
two between each other. Model 3 is qualitatively the same, but the 
amplitude is twice as big and the oscillations one third less frequent 
with respect to model 1. These models allow a comparative study 
on the influence of the amplitude and the frequency of the oscilla- 
tions. Models 4 and 5, despite being described by the same func- 
tional form, have a very long period and oscillations are not even 
visible throughout cosmic history. Models 3 and 5 are also charac- 
terized by a crossing of the phantom barrier, w — —1, a feature 
that is marginally find to best fit the luminosity distance of SNe la. 
Finally, model 6 shows tiny oscillations only at recent times, while 
for z > 10 the function w{z) approaches a constant. By comparing 
model 6 with the others we can draw conclusions about the impor- 
tance of the oscillations at early times. 

3.1 Homogeneous background analysis 

We now explore in detail the redshift evolution of functions re- 
lated to the homogeneous background for the various oscillating 
quintessence cosmologies presented above. In the upper panel of 
Figure |2] we show the ratio between the matter density parameters 
in the six dynamical dark energy models considered here and the 
same function in the concordance ACDM cosmology, as a func- 
tion of redshift. In the lower panel of the same Figure we display 



-0.2 
-0.4 




1 1 1 . 1 1 


Model 1 
Model 2 
Model 3 
Model 4 
Model 5 
Model 6 


1 . II 
- 


-0.6 










-0.8 
-1 


\ •'••<./ ' h \ 


..JS'". 






-1.2 


. . 1 . 


.1 , . . 1 


. . 1 . 


1 . 



10° 10^ 10^ 10^ lO'* 10^ 10^ 



1+z 

Figure 1. The redshift evolution of the equation of state parameter for the 
oscillating dark energy cosmologies analyzed in this work. The red dashed, 
blue short-dashed and cyan dotted curves show the models 1 , 2 and 3 respec- 
tively. The orange dot-dashed curve shows model 4, while the dark-green 
dot-short-dashed and the light-green dot-dotted lines represent models 5 and 
6, respectively (see Table[T]for more details). 

the corresponding redshift evolution for the dark energy density 
parameters. As one could naively expect, the amount of matter at 
early times is the same in all models, consequence of the fact that at 
high redshift the dark energy contribution becomes negligible, and 
hence the Universe always behaves as an Einstein-de Sitter (EdS) 
cosmology. Differences become significant at z < 5 and are at most 
at the level of ~ 10 — 15%. It is worth noticing that in no circum- 
stances we see full oscillations in the density parameters, implying 
that one integral over the cosmic history is enough to smooth out 
most features of wlz). 

It is worth noting that model 6 and, especially, model 4 dis- 
play very little difference with respect to the concordance case, an 
instance that will show up time and again throughout the discus- 
sion of our results. The other models instead, with the exception 
of model 5 that sees a reduction in the abundance of dark energy, 
show a substantially higher amount of dark energy at early times 
than the cosmological constant case. At 2 ~ 20 the difference in 
the dark energy density parameter is of one order of magnitude or 
more. It should be recalled however that at such high redshifts the 
contribution of dark energy to the expansion history of the Universe 
is negligible anyway. 

A different perspective on the same results is given by exam- 
ining the Hubble parameters (that are, the expansion rates) for the 
various cosmologies, shown in the upper panel of Figure[3] For flat 
universes, the Hubble parameter can be written as 

H{a) = HoE{a) = ^^o^^ + t^q.o g{a) , (II) 

where 17(a) is defined as 

Sr(a) = exp ( -Sy -j^^da J , (12) 

and Q,q is the dark energy density. The behavior of the expansion 
rates is very similar to that of the matter density parameters, which 
is expected because dark energy comes to dominate the evolution 
of the Hubble function only at very low redshift, where differences 
between different models tend to vanish. 

Since oscillations are not present in the redshift evolution of 
the matter and dark energy density parameters, the same holds true 
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Figure 2. The redshift evolution of the density parameters. In the upper 
(lower) panel we show the ratio of the matter (dark energy) density parame- 
ter in the six oscillating dark energy cosmologies studied in this work to the 
corresponding function in the fiducial ACDM model. Line types and colors 
are as in Figure [T] 



for the expansion rate in our quintessence cosmologies. However, 
considering the derivative of the expansion rate with respect to 
the scale factor (dubbed deceleration parameter, see Dunajski & 
Gibbons 2008 and references therein), shown in the lower panel 
of Figure [S] we observe some partial indication of oscillations, in 
that the oscillating dark energy deceleration parameter crosses the 
ACDM one at least once for models 1, 2, and 5. However, since 
the overall pattern looks the same for all models, it is likely not di- 
rectly connected with dark energy oscillations. Moreover, oscillat- 
ing quintessence models introduce absolute differences in the de- 
celeration parameter of at least ~ 20 — 30% with respect to the 
fiducial case, which is a quite significant effect. 

In the upper (lower) panel of Figure |4] we present the ratio of 
the luminosity distance (age of the Universe) in the various dynam- 
ical dark energy cosmologies to the same quantity in the ACDM 
model, as a function of redshift. In both cases we see that differ- 
ences are at most of the order of ~ 4% and predominantly located 
at relatively low redshifts (z < 2), although the luminosity distance 
shows a ~ 2% deviation even at arbitrarily high redshifts. This 
fact is expected, since both the luminosity distance and the age of 
the Universe are suitable integrals over some function of the Hub- 
ble parameter, which also shows most differences at low redshift. 
Comparing the upper panel of Figure [4] with the luminosity dis- 
tances inferred by SN la Union2 data l lAmanuUah et al.|[2010l) we 
see that differences induced by the oscillating dark energy models 
at low redshifts are at the same level of the systematic errors in the 
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Figure 3. Upper (lower) panel. Ratio of the Hubble parameter (derivative of 
the Hubble parameter with respect to the scale factor) in the six oscillating 
dark energy models considered in the present paper to the same function in 
the fiducial ACDM cosmology. Line styles and colors are as in Figure[T] 



measurements as well as of the intrinsic scatter around the best fit. 
This, together with the very slight deviations in various cosmolog- 
ical functions shown in previous Figures lead us to conclude that 
the oscillating quintessence cosmologies studied in this work are 
not distinguishable from the concordance model by current geo- 
metrical probes. This is perfectly consistent with previous works, 
since the parameter values that we adopted are indeed chosen so as 
to reproduce some particular geometrical tests. 



3.2 Redshift drift 

An important cosmological test related with the expansion history 
that has not been considered in the past, but might reveal itself valu- 
able in the near future is the so called redshift drift, that represents 
the variation of the cosmological redshift of a source du e to the 
expansion of the Univer se (Balbi & Ouercellini 2007; Lis ke et al. 
'20081; lUzan et alj |2008| ; Uain & JhingarJ 12010 ; Arailjo & Stoege] 
2010). Let us indicate with is the time of emission of an electro- 
magnetic signal from a source, and with to the time of observation 
of the same signal. The cosmological redshift of the source is then 
defined as 



1 + Zs 



a{U) 



(13) 



After a time interval Afo has passed for the observer, correspond- 
ing to an interval Ais for the source, the change in the source red- 



6 F. Pace et al. 



1.12 
1.1 
1.08 
1.06 
1.04 
1.02 
1 

0.98 
0.96 

1.05 
1.04 
1.03 
1.02 
1.01 
1 

0.99 
0.98 
0.97 
0.96 



Model 1 
Model 2 
Model 3 
Model 4 
Model 5 
Model 6 



10 

z 



15 



20 



1 1 


Model 1 

Model 2 - 

Model 3 
Model 4 

Model 5 

Model 6 










1 1 


1 



10 

z 



15 



20 



Figure 4. Upper (lower) panel. Ratio of the luminosity distances (ages of 
the Universe) in the oscillating dark energy models considered here to the 
same quantity for the tiducial ACDM cosmology, as a function of redshift. 
Line styles and colors are as in Figure[T] 



shift can be estimated by expanding at first order the previous equa- 
tion, 



a(io) - a{ts) 



(14) 



By substituting H[a) = d/awe obtain the expression for the red- 
shift drift 



Azs ~ HoAto [1 + 2s - E{zs)] , 



(15) 



where E{zs) = H{zs)/Ho. 

Using the variation in the cosmological redshift, it is also pos- 
sible to determine the variation in the recession velocity of the 
source. 



cAZa 
1 + 2a 



We can therefore write 



E{zs 



(16) 



(17) 



Since all the cosmological properties of the model at hand are en- 
coded into the Hubble function, we see that we can use the time 
variation of the redshift in order to reconstruct the expansion his- 
tory of the Universe. 




Figure 5. The difference between the redshift drifts in the six oscillating 
quintessence cosmologies considered in this work and the some quantity in 
the fiducial ACDM model, as a function of redshift. Line styles and colors 
are the same as in Figure[T] 



In Figure[5]we show the difference between the redshift drifts 
in the oscillating dark energy models and the same function in the 
concordance ACDM cosmology, per unit of observed time and nor- 
malized by the Hubble constant. We decided to plot the differences 
instead of ratios in this case in order to avoid divergences, due to 
the fact that Azs goes to zero when 1 + Zs = E{zs). Similarly 
to previous Figures, also in this case models 4 and 6 behave very 
similarly to the ACDM cosmology. The redshift drift is systemat- 
ically higher for the model 5 while it is systematically lower for 
the first three oscillating cosmologies. At high redshifts the differ- 
ences between the models decrease since all the Hubble parameters 
converge to the EdS behavior. Since all the cosmological informa- 
tion is encoded in the Hubble expansion function, no oscillations 
appear in this case as well, although a slight wiggle is visible for 
model 1 at 2 ~ 2. As for the pe rspective of realist ically measur- 
ing the redshift drift, according to lLiske et al.l ( l2008h . peculiar mo- 
tions are negligible (~ 10~^ cm s~^) and with a temporal base- 
line of 20 years it will be possible to determine the existence of 
the cosmological constant at 3.1 g observing distant quasars. Fol- 
lowing Balbi & Ouercellinj | |2007|) we also notice that variations 
in the recession velocity of the sources are bigger than the error 
bars forecasted by Montecarlo simulations, therefore with a suffi- 
ciently long baseline, it should be possible to discriminate between 
oscillating quintessence and cosmological constant, at least for the 
models showing more significant deviations. 



4 RESULTS 

In this Section we present results concerning the structure forma- 
tion in the oscillating quintessence cosmologies described in Sec- 
tion[3] We studied several aspects of structure formation and in par- 
ticular we focused our attention on the growth factor, the linear and 
non-linear overdensities derived from the spherical collapse model, 
the mass function of cosmic structures, the power spectrum of cos- 
mic shear and the ISW effect. We now proceed to describe each one 
of these observables in detail. 



4.1 Growth factor 

In Figure |6] we show the growth factor normalized by the scale 
factor, D+[a)/a, as a function of redshift for the six oscillating 
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dark energy cosmologies described above, plus the fiducial ACDM 
model. Cosmological observables sensitive to the growth factor 
include cosmic shear, ISW effect and the Rees-Sciama effect, all 
discussed later on. In Figure |6] the growth factor is normalized to 
unity at 2 = 0. As can be seen, differences between the oscillat- 
ing quintessence models and the ACDM cosmology (solid black 
curve) are at most of ~ 10% (for model 3), while, in agreement 
with previous Figures, model 4 does not show any appreciable dif- 
ference from the concordance model. It is interesting to note that, 
while w{z) for model 4 shows indeed very little variation up to the 
last scattering redshift due to the very large period of its oscillations 
(see Figure[TJ, model 6 shows even smaller time evolution, yet its 
effects on cosmological functions are larger. This implies that high 
frequency oscillations 'mw[z), albeit with a very small amplitude 
and limited time extent have more of an effect on the expansion his- 
tory (and structure formation too, see later Subsections) than larger 
oscillations with a low frequency. This is because low-frequency 
oscillations cancel integral contributions more effectively. 

Figure[6]also shows that for models 5 and 6 the growth factor 
is smaller than for the cosmological constant case. This can be un- 
derstood by the following argument. From the lower panel of Fig- 
ure |2] we see that for these two models the amount of dark energy 
is smaller than for the cosmological constant case at all redshifts. 
This means that the Hubble drag is less effective in the former mod- 
els, and hence the growth of structures (at least at the linear stage) 
is easier. Since the amplitude of density fluctuations at 2 = is 
almost the same amongst all the dark energy models considered 
here (it differs by a factor proportional to the critical overdensity 
for spherical collapse, that however is only slightly changed in the 
case of dynamical dark energy, see below), the growth factor must 
be smaller in order to match the amplitude of fluctuations at early 
times. 

The growth factors depicted in Figure |6] do not present any 
sign of oscillations, not even if we consider their ratios with respect 
to the concordance cosmology case. By rewriting Eq. ^ using the 
scale factor instead of cosmic time as the independent variable, we 
obtain 



5" + 



3 



S' 



3 n„ 



(18) 



As shown above, the derivative of the Hubble function presents 
mild signs of oscillations, however since this is a second order dif- 
ferential equation, the solution involves a double integral over the 
scale factor, that efficiently smoothes out any fluctuation in the co- 
efficients. 

Finally, we also estimated the logarithmic derivative of the 
growth factor with respect to the scale factor, f{Q,n-i{a)) = 
d In (a) /d In a. It has been shown that in a broad range of cos- 
mologies f{yim[a)) ~ ^Ziio)^ sri empirical relation that we re- 
trieve for the oscillating dark energy cosmologies as well. Devia- 
tions with respect to the ACDM expectation however are smaller 
than for the growth factor itself. 

4.2 ISW and Rees-Sciama effects 

The Sachs-Wolfe (SW) effect ( ISachs & Woif^ligeTl) describes the 
effect of gravitational potentials on the CMB anisotropy at the last- 
scattering surface. Photons travelling to an observer encounter vari- 
ations in the gravitational potential caused by variations in the mat- 
ter density. Photons climbing out a potential well will be gravita- 
tionally redshifted and this will make the region they come from ap- 
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Figure 6. The growth factor as a function of the redshift. The black solid 
line represents the concordance ACDM cosmology, the red dashed, blue 
short-dashed and cyan dotted curves show the models 1, 2 and 3 respec- 
tively. The orange dot-dashed curve shows model 4, while the dark-green 
dot-short-dashed and the light-green dot-dotted represent models 5 and 6, 
respectively. 
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Figure 7. Redshift evolution for the quantities describing the ISW (up- 
per panel) and the Rees-Sciama (lower panel) effects. For the ISW effect 
we plot the function d(D^(a)/a) /da, while for the RS effect we plot 
d {Dl{a)/a) /da. In both cases, the result for the fiducial model is sub- 
tracted from the corresponding quantities for each of the six oscillating dark 
energy cosmologies. Line types and colors are the same as in Figure [T] 
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pear colder. Together with this gravitational effect, one has to take 
into account the time-dilation effect induced by the perturbations: 
we see the photons as coming from a different spatial hypersurface 
(labeled by a different scale factor a{t)). 

The Integrated Sachs-Wolfe (ISW) effect is based on the same 
principle, only it is given by the gravitational redshift occurring 
as photons travel through the large scale structure to reach an ob- 
server at present time. The ISW effect arises only recently in the 
cosmic history, as dark energy starts dominating the expansion of 
the Universe. This means that a non-vanishing ISW effect indicates 
by itself presence of dark energy if the model is spatially flat, as 
indeed it is assumed in our case. The ISW effect is sensitive to the 
derivative of the growth factor, d {D+ (a) /a) /da that vanishes for 
an EdS universe where D+(a) = a . It was detected for the first 
time bv lBoughn & Crittenden j2004) using X-ray c luster catalogs. 
The Rees-Sciama (RS) effect ( Rees & Sciamall968h is very similar 
to the ISW effect, only it refers to the gravitational redshift induced 
by non-linear structures only, and hence it is active on much smaller 
scales. It is mainly sensitive to the function d (_D+ (a) / a) /da. 

In Figure |7] we show the difference of the functions probed by 
the ISW effect (upper panel) and by the RS effect (lower panel) 
for the six oscillating dark energy models studied here to the same 
quantities evaluated in the reference ACDM cosmology, as a func- 
tion of redshift. As can be seen, the ISW effect is preferentially 
modified at low redshifts, either positively or negatively, except for 
models 4 where no differences from the ACDM model are seen. At 
high redshifts deviations with respect to the cosmological constant 
case tend to disappear, since all the models are very well approxi- 
mated by an EdS universe. The models showing the largest effect 
are models 3 and 5, which are the ones having the largest differ- 
ences in the growth factor. On the other hand, hints of an oscil- 
latory behavior with redshift are seen for model 1, which has the 
highest frequency in w{z) among those considered here. All in all, 
since differences between different models can be quite substantial, 
high redshift observations could be used in principle to discrim- 
inate oscillating quintessence cosmologies. Examples come from 
cross-correlating galaxies, radio so urces or hard X-ray so urces and 
CMB temperature fluctuations, (see Fosalba et al . 2003; Nolta et al.l 
2004 IScranton et al]|2003l : iBoughn & CrittendenI 12004 ; lAfshordil 
2004) . However, one should keep in mind that at high redshift, 
where differences are more marked, the ISW effect itself tends to 
disappear. 

As for the RS effect, deviations with respect to the cosmolog- 
ical constant expectations are of the same order of magnitude as 
the ISW effect and differences do not vanish at high redshifts, but 
reach a somewhat constant value, the exact one depending on the 
specific model. There is however a very specific redshift pattern ac- 
cording to which the difference with respect to the fiducial ACDM 
cosmology gets reversed at z ~ 2 (except for model 4, that is basi- 
cally identical to the concordance cosmology). Hence, combining 
RS effect observations at low and high redshift can improve the 
discrimination between the models. 



4.3 The characteristic overdensities 5c and Av 

In this Subsection we present results regarding the time evolution 
of the linear density contrast parameter for spherical collapse Sc 
and of the virial overdensity Av. The main results are reported in 
Figure[8] In the upper panel we show the time evolution of Sc while 
in the lower panel we present the time evolution of Av. Line types 
and colors are the same as in Figure (6] to which we refer for a 
detailed explanation. The first thing to note is that, contrary to ex- 
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Figure 8. The redshift evolution of the linear density contrast parameter 
<5c (upper panel) and of the virial overdensity parameter Av (lower panel) 
for the six dynamical dark energy models and for the reference ACDM 
cosmology. Line types and colors for the different models are as in Figure[6] 



pectations, the function Sc{z) does not perfectly converge to the 
EdS value of 5c — 1.686 at high redshift. This is a problem of nu- 
merical convergence related to the oscillatory nature of w{z), and 
that is better explored in Appendix |A] This fact obviously implies 
that we should not use this 5c{z) at 2 > 5 — 6. However, all the 
cosmological tests that we propose in the following that make use 
of this function are limited to substantially lower redshifts, hence 
they should be unaffected by this issue. 

As mentioned above, differences between the 5c{z) computed 
in different cosmologies are very mild, being at most of ~ 1% at 
z < 2. It is apparently a generic feature of cosmological models 
displaying a dynamical evolution of the dark energy density that 
the spherical collapse parameters ar e only slightly modified with 
respect to the fiducial ACDM case jPace et al. 2010). The lower 
panel of Figure [8] shows results for the virial overdensity param- 
eter Av(z). In order to evaluate it we used the prescription of 
IWang & Steinhardtl ( Il998h . In this case for z > 6 — 8 all the mod- 
els behave almost exactly as the ACDM cosmology. Differences 
between different models are of the order of a few percent and are 
mostly evident at 2 < 3. We tr ied to evaluat e Av(2:) adopting a dif- 
ferent prescription, detailed in lWand ( |2006|) . As it turns out, differ- 
ences between different models are very similar to those obtained 
by using other recipes. 
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4.4 Mass function 

An observable quantity depending crucially on the growth factor 
D+(z) and on the linear overdensity threshold for collapse 5c(z) 
is the mass function of cosmic structures n{M, z), representing the 
number of dark matter halos per unit mass and per unit comoving 
volume present at a certain redshift. Integrals of the mass function 
over mass can be observed directly by using large cosmological sur- 
veys, provided their selection function is well understood. Specifi- 
cally, the cumulative number density of cosmic structures above a 
certain limiting mass Mmin (that will depend on the specific sur- 
vey) at redshift z is simply given by 

N{>M^in,z)= dMn{M,z). (19) 

The mass variance is another key ingredient for the mass func- 
tion formalism, and is identified by 

<TM = 2^ k''T^{k)WUk)Po{k)dk . (20) 

In Eq. J20b Po (k) represents the primordial matter power spectrum, 
T{k) is the matter transfer function, while WR{k) is the Fourier 
transform of the real space top-hat window function. Since the only 
difference induced by oscillating dark energy in the primordial mat- 
ter power spectrum is given by the different normalization erg, that 
as we shown above is very slight, we expect only minor differences 
in the mass variance as well. In this work, to evaluate the mass 
funct ion, we used the Sheth-Tormen expression ( Sheth & Tormei] 
[19991) . 

In Figure |9] we show the cumulative mass function for the six 
oscillating quintessence cosmologies considered in this work at dif- 
ferent redshifts, divided by the corresponding quantity evaluated 
for the reference ACDM model. From this Figure we note first of 
all that there are no differences between different models at 2: = 0. 
This is due to our choice of the normalization, namely that the ra- 
tio Sc{z — 0)/(7g should be conserved for all models. Significant 
differences start instead to appear at higher redshifts, where we no- 
tice that the six models can be broadly divided into two groups. 
Models 1, 2, and 3 show more objects with respect to the ACDM 
case, while models 4, 5, and 6 show less objects. These differences 
are consistent with the latter group having a lower growth factor 
with respect to the former and a higher critical linear overdensity 
Sc{z), which makes more difficult for density perturbations to col- 
lapse into bound structures. As one could naively expect, deviations 
increase with increasing redshift and are most notable in the very 
high mass tail, since rare events are very sensitive even to small 
fluctuations in the expansion history. As an example, differences in 
the abundance of massive galaxy clusters A/ > 5 x lO^*M0//i 
range from ~ 5 — 10% at z = 0.5 up to ~ 30% at 2 = 1. 
Large galaxy groups that can be expected to be found at 2 = 2 
(M > 5 X IO^^Mq/Zi) are ~ 20 - 30% more abundant in models 
belonging to the first group, and up to ~ 40 — 50% less frequent 
in models of the second group. We also note that models showing 
an enhancement in the abundance of cosmic structures are those 
having more dark energy at high redshift. 

It should be noted that up to now we only considered the co- 
moving number density of objects, that is we did not take into ac- 
count possible effects deriving from changes in the cosmic volume 
induced by the presence of oscillating dark energy. We will include 
this additional ingredient shortly. For the time being, we can con- 



clude that the impact of oscillating quintessence on the counts of 
cosmic structure can be quite substantial, especially at high red- 
shift, thus implying that a detection might be possible with future 
large cluster surveys. 

In order to establish a more direct link with observations, we 
forecast the redshift distribution of galaxy clusters that, in each of 
the cosmological models considered in this work, will be observed 
by upcoming wide field surveys with cluster selection based both 
on Xray and weak-lensing data. In order to do that we need to define 
a redshift dependent minimum mass for the observed objects, and 
integrate the mass function above that threshold. We assume that 
precise estimates for the masses of these objects will be available, 
which is realistically expected if a robust multiwavelength follow- 
up will be performed. 

The first survey that we consider is a wide field X-ray sur- 
vey on the model of the upcoming eROSIT^H one. In order to 
determine the minimum mass of clusters that will compose the 
X-ray catalog we need a scaling relation between the observable 
at hand (in this case the X-ray flux) and the true mass. First of 
all, knowing the redshift of the cluster, the measured bolomet- 
ric flux can be related to the intrinsic bolometric luminosity as 
I/(M, 2) = 47r/(A/, z)d'^{z), where di,{z) is the luminosity dis- 
tance (see upper panel of FigureQ. In reality the X-ray bolometric 
flux is almost never measured, rather the X-ray photon counts in a 
certain energy band are used. For the specific case of an eROSITA- 
like X-ray survey, we adopted the band [0.5, 2.0] keV, where the 
threshold flux is expected to be /min ~ 3.3x 10~^^ erg/(s cm^). In 
order to estimate the band flux for a cluster of a given mass placed 
at a given redshift we modeled the related intra-cluster m edium us- 
ing a Raymond-Smith plasma model ( Raymond & S mi thil 19771) as 
implemented in the Xspec softw are package ([A maud 1996), with 
metal abundance Z = 0.3^0 jFukazawa et al.. 1998; Schi ndle3 
1 19991) . The plasma model has been normalized so as to repro- 
duce the bolom etric luminosity expected from the s caling relation 
adopted bv Fedeli. Moscardini & Bartelmanj ( l2009f) , namely 

L{M, 2) = 1.097 X lO**^ erg/s {ME{z)Y-^^'^ , (21) 

where the mass must be inserted in units of IO^^A/q/Zi. This rela- 
tion results from the combination of two scaling laws, one relating 
the mass to the X-ray temperatu re, and the other relating the tem- 
perature to the luminosity. See iFedeli. Moscardini & BartelmannI 
(2009) for additional details. 

The second case that we consider is representative of clus- 
ter catalogs selected through their weak-lensing signal. Massive 
galaxy clusters can be selected as high S/N peaks in the weak lens- 
ing map produced by weak-lensing surveys. The S/N strength will 
also be used as a proxy for their mass, although a robust multi- 
wavelength follow-up program will be necessary in order to have 
more precise estimates. For the minimu m cluster mass entering this 
catalog, we adopted the calculations of iBerge. Amara & Refregien 
In particular, we refer to their Figure 1 , wh ere they present 
the selection function for a Eiiclid-liks surve >0 dLaureiijlioO^: 
iLaureiis et al.l20llb in the mass-redshift plane, assuming a number 
density of background galaxies of fig = 40 arcmin"^. We consid- 
ered the contour referring to a S/N threshold of 5, threshold that 
was shown to be a good choice to minimize spurious detections in 
the weak lensing maps (see iPace et aLj2007i) . 



^ http://www.mpe.mpg.de/erositay 

^ http://sci.esa. int/science-e/www/area/index.cfm?fareaid=102 
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Figure 9. Cumulative comoving number density of objects with mass exceeding M at different redshifts. Ratios with respect to ACDM expectations are 
shown. Selected redshifts are z = Q (upper left panel), z = 0.5 (upper right panel), z = 1 (botton left panel) and z = 2 (bottom right panel). Line styles and 
colors for the different models are as in Figure[T] The upper left panel is unity by normalization, as explained in the text. 



In Figure [To] we present our results for the minimum mass of 
the catalogs. In the upper panel we show the minimum mass for the 
adopted weak-lensing survey (sohd line) and X-ray survey (dashed 
line) for the ACDM model. We see that as expected both mini- 
mum masses increase with redshift in order to have the same flux 
or S/N ratio. The X-ray mass increases from few times 10^'^ Mq /h 
at 2 ~ till ~ 6 X 1O"M0 /ft at z = 2. As evident, since the flux 
limit is constant in redshift, the redshift dependence of the mass can 
be very well approximated by a parabola, therefore compensating 
the increase of the luminosity distance (entering quadratically in 
the relation between flux and luminosity). The minimum mass for 
a weak-lensing survey is increasing much faster with redshift since 
the lensing efficiency drops very fast to zero if the lens approaches 
the sources. 

In the lower panel we show the minimum cluster mass for the X-ray 
catalog in each dark energy cosmology considered here, divided by 
the same quantity estimated for the reference ACDM model (lower 
panel). We see that differences in the minimum observed mass are 
strongly related to the differences in the Hubble function and in the 
luminosity distance, as one might naively have expected. Specifi- 
cally, since the minimum mass depends on the square of the lumi- 
nosity distance, even small variations of the latter turn out to affect 
the former at the level of ~ 10%. 

In Figure[TT]we present the redshift distributions expected for 
the weak-lensing survey (upper panel) and the X-ray survey (lower 
panel) we considered, respectively. The redshift distribution is de- 
fined as 



dV{z) 
dz 



N{> Mmin{z),z), 



(22) 



where Mminiz) is the minimum observed mass for the survey at 
hand and dV{z)/dz is the comoving volume element contained 
in the unit redshift. The first thing to note is that the deviations 
in the redshift distributions induced by the presence of oscillating 
dark energy are quite substantially different for an X-ray survey and 
and a cosmic shear survey. This is likely due to the fact that these 
surveys have remarkably different selection functions, that sample 
quite distinct regions of the mass-redshift plane, where the impact 
of oscillating dark energy is necessarily different. Let us consider 
first the X-ray eROSITA-like survey. In this case the redshift distri- 
butions computed for different cosmologies are almost identical at 
very low redshift, while substantial deviations are visible at higher 
redshift. Specifically, the abundance of clusters is incremented by 
up to ~ 20% at z ~ 0.8 for models 1, 2, and 3, and decreased by 
the same amount or more for the other models. 

For the weak lensing survey the situation is totally different. 
Cosmological models differ from each other by up to ~ 5 — 10% 
already at relatively low redshift. Deviations from the ACDM cu- 
riously vanish at 2 ~ 0.6 for all models, and then grow again, but 
with opposite sign, for higher redshifts. In addition to the differ- 
ent selection function, one additional difference between the X-ray 
survey and the cosmic shear survey considered here is that in the 
former case the minimum mass depend on cosmology, while in the 
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Figure 10. emphUpper panel. Minimum mass for weak-lensing survey 
(solid line) and X-ray survey (dahsed line) for the reference ACDM model. 
Lower panel. Minimum mass for the eROSITA cluster catalog as a function 
of redshift, presented as ratios between the oscillating dark energy models 
and the ACDM cosmology estimates. 



latter case it is model independent since what we measure is di- 
rectly related to the mass of the cluster. 



4.5 Dark matter power spectrum 

An important tool that can be used to infer the statistical proper- 
ties of a cosmological model is the fully non-linear matter power 
spectrum. Observationally, this can be estimated both by using 
the distribution of tracers such as galaxies and galaxy clusters 
(provided their non-linear bias is understood) and through cosmic 
shear (that however returns only a projected version of the three- 
dimensional spectrum). The matter power spectrum can be stud- 
ied in the non-linear regime by means of numerical N-body sim- 
ulations or semi-analytic prescriptions that are fitted against them 
jPeacock & Doddslll996l : Ismith et alj|2003h . The shortcoming of 
such fitting formulas stays in the fact that they have limited tests 
of validity, usually restricted to the concordance ACDM comology 
and scales fc < 50 — 100 h/Mpc at z = 0. 

An alternative approa ch, that we exploi ted, is based on the 
halo model developed by iMa & Frvl ( |2()0(]|) and ISeliakI ( l200(]h 
(see ICoorav & Shetl]|2002 ! for a comprehensive review). The halo 
model has a physical motivation, and relies on ingredients, such 
as the average dark matter halo density profile and mass function, 
whose universality is much better established than for the non- 
linear matter power spectrum. Additionally, using this formalism 
the matter power spectrum calculation can be pushed to very small 
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Figure 11. Cluster redshift distribution, as a function of redshift, presented 
as ratios between the oscillating dark energy models and the ACDM cos- 
mology estimates. Upper panel. Number of objects above the minimum 
mass for a Euclid-\i\ie survey. Lower panel. Number of objects above the 
minimum mass for an eROSITA-like survey. Line styles and colors for the 
different models are as in Figure[T] 



scales, in principle as small as the structure of cold dark matter ha- 
los has been studied. Within the halo model the full power spectrum 
is given by the sum of two terms: one (the 2— halo term) dominates 
on large scales and it depends on the correlation of individual halos; 
the other (the 1— halo term) dominates on small scales and it is sen- 
sitive to the inner structure of the halos. The ingredients needed in 
order to apply the halo model formalism are the halo mass function, 
the halo linear bias, and the halo mass density profile. 

The behavior of the mass function in the various oscillating 
dark energy cosmologies considered here has been described in 
the previous Subsection. For the average internal structure of dark 
matter halos we adopt the Navarro. Frenk & White, C1997.) (NFW 
henceforth) density profile, 



p(r) = 



(r/r,)(l + r/r-e)2 



(23) 



where ps is a density scale while rs is the radius at which the log- 
arithmic slope of the profile equals —2. The parameter is re- 
lated to the virial radius i?v of the structure by the concentration, 
c = Rv/ra- The concentration is in turn depending on the total 
mass M of the object, in that dynamically younger structures have 
on average larger masses and lower concentrations. We adopt the 
following concentration-mass relation. 
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c(M, z) 



10 

l + z 



M 



M,{z) 



(24) 



inspired bv lHuffenberger & SeliakI l l2003h and tested to give good 
agreement with fits to N-body simulations in the framework of the 
ACDM cosmology. In the previous equation A/, [z) represents the 
characteristic collapsing mass at a given redshift, defined implicitly 
as D+{z)(JM^ = 5c{z). 

As for the large-scale bias, which is needed in the 2— halo 
part of the power spectr um, we adopt the prescription by 
ISheth. Mo & TormenI 1 I2OO1I) , 



b{M,z) = l + a 



Sc{z 



+ 



_2p_ 



Dl{z)ol, 5^{z) 

{D+{z)aMf^ 



{D+{z)gmYp + {•/^5,{z)Yp 



(25) 



where a = 0.75 and p = 0.3. 

In Figure [12] we show the ratio of the matter power spectrum 
computed in each of the six oscillating dark energy models ex- 
plored here to the same function in the concordance ACDM cos- 
mology, as a function of wavenumber. In the upper panel we plot 
the ratio at 2 = while in the lower panel we consider the ratio at 
z = 0.5. It is immediately evident that at the present time all cos- 
mologies tend to have the same non-linear matter power spectrum 
at small scales. This is due to the fact that in the concentration-mass 
relation that we adopted (Eq. I24l l all the cosmology dependence is 
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Figure 13. The effective convergence power spectrum for each of the six 
oscillating dark energy cosmologies considered in this work. We show ra- 
tios with respect to the ACDM power spectrum. Line styles and colors for 
the different models are as in Figure [T] 



encapsulated in the characteristics non-linear mass M,{z). How- 
ever, due to the fact that D+{z = 0) = 1, and Sc{z = 0)/as 
is the same for all cosmological models, the concentration for a 
given mass at z = is also unchanged. Hence, since dark matter 
halos share always the same inner structure, they produce an iden- 
tical dark matter correlation function at small scales. The situation 
changes at 2 > 0. In fact, as we can see from the lower panel of Fig- 
ure [12] the matter power spectra now differ from the ACDM case 
at all scales. Differences are however at the level of ~ 1% at most 
and model 4 in particular does not show strong differences from the 
reference model. One might argue that the inner structure of dark 
matter halos should change at least a bit due to oscillatory dark en- 
ergy, even at 2 = 0. Since however we have no indication on how 
this is expected to happen, we stick to our original choice. At very 
large scales the power spectrum approaches the linear primordial 
one, hence the deviations with respect to the fiducial cosmology 
are induced only by differences in the normalization as (squared), 
that are always below 1%. 



4.6 Cosmic shear 

As mentioned above, the power spectrum of cosmological weak 
lensing is a projected (and weighted) version of the fully non-linear 
matter power spectrum. What is actually measured is the effective 
convergence power spectrum (equivalent to the shear and reduced 
shear power spectra), given by 



iKixY 



(26) 



where a(x) is the scale factor, x the comoving distance up to scale 
factor a, and /if (x) the comoving-angular diameter distance which 
depends on K, the spatial curvature of the Universe. The integral in 
the previous equation formally extends up to the horizon size xh, 
however since the number density of sources (see below) drops to 
zero much before that, the integral can be effectively truncated at 
2 ~ 10. 

The function W{x) is an integration kernel depending on 
n{z{x)), the redshift distribution of background sources. The ker- 
nel can be written as 
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dx' 



(27) 



The functional form for the redshift distribution of the sources has 
to be inferred with the help of observations. In the following wede- 
cided to assume the source distribution derived bv lFu et aU i2008h 
using data of the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS). The distribution takes the functional form 



n{z) = A 



z'> + c 



(28) 



with parameter values a = 0.612, b - 
normalization constant A is given by 



.125, and c = 0.62. The 



A' 



z" + C 



dz . 



(29) 



We checked our results against variations of the source red- 
sh ift distributio n, verifyi ng that adopti ng instead the distribution 
of Brainerd, Bl andford & Smaill ( Il996h (see also lEfstathiou et al.l 
ll99lUSmail&'Dickinsoijll995l) or lseniamin et al.(j2007l) changed 
very little the subsequent cosmic shear results. 

In Figure [T3] we show the effective convergence power spec- 
tra for the various oscillating dark energy models explored in this 
work, divided by the same quantity estimated in the framework of 
the concordance cosmology. We can observe that models 4 (yellow 
dot-dashed curve) and 6 (cyan dot-dotted curve) have negligible 
differences with respect to the ACDM model, while for the other 
models deviations can reach up to ~ 10 — 15% at very small an- 
gular scales. The fact that the largest deviations are visible at small 
scales highlights how the matter power spectrum is indeed affected 
by oscillating quintessence at those scales if one considers z > Q. 

Let us now make the case for a possible detection of oscillating 
dark energy more specific, by considering the S/N ratio for such a 
detection at a fixed multipole. This can be written as 



N 



dACDM 



(30) 



where AP^^^^{£) is the Gaussian statistical error on the power 
spe ctrum in the framewo r k of the conco r dance co s molog y. Follow- 
ing iKaiseJ ( 1 1991 Il998h : ISeliakI ( Il998l) ; iHutereil ( 120021) , the latter 
can be evaluated as 
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Figure 14. The S/N ratio for discriminating between the concordance cos- 
mology and each one of the six dynamical dark energy models considered in 
this work, as a function of multipole. Line styles and colors for the different 
models are as in Figure[T] 



in our analysis, kick in, and baryonic physics cannot be neglected 
anymore as well. Given this, in Figure [T4l we show the S/N ratio 
for the effective con v ergen ce power spectrum. As also found by 
iFedeli & Moscardinil ( l201oh when studying the impact of primor- 
dial non-Gaussianity on the weak lensing power spectrum, we ob- 
serve that at intermediate scales, 100 < £ < 1000, S/N ~ 0.5 - 1 
for four out of the six models we studied, while for very low or 
very high multipoles the S/N ratio tends to vanish. This means that 
for these models it would be sufficient to sum the S/N ratio over 
a relatively limited number of intermediate multipoles in order to 
have a significant detection of oscillating dark energy over the con- 
cordance cosmology. In models 4 and 6 however, the S/N ratio is at 
most at the level of ~ 10^'^, hence these models would be very dif- 
ficult to distinguish from the ACDM cosmology by using a cosmic 
shear survey. 

A maneuver usually capable of increasing the discriminatory 
power of cosmic shear consists in subdividing the source redshift 
distribution in a certain number of bins (usually 5 is the maximum 
number that gives appreciable improvement, see ISun et al1l2009h . 
computing the effective convergence power spectrum given by each 
bin, and then combining the vario us spectra together. This expedi - 
ent is dubbed lensing tomography jHil 19991 ; iTakada & Jainll2004h . 
Specifically, the cross spectrum between the two source redshift 
bins i and j is just a straightforward generalization of the power 
spectrum defined above in Eq. l |26l l. 



(2f + l)A^/,ky 



(31) 



where fig is the average surface number density of observed galax- 
ies, /sky is the fraction of sky area surveyed, and 7 represents 
the rms intrinsic shape noise for the average galaxy. By assum- 
ing the specifications of a future typical weak lensing surve y we 
set n = 40 arcmin"'^, /sky = 1/2 and 7 — 0.22 (see lZhang et al.l 
I2OO9I) . To specif y A^, the binwidth over whic h the spectrum i s aver - 
aged, we follow lTakada & Bridld i2007h and lTakada & Jaiij ^20091) 
and adopt the value A^ = 1. A surv ey like Euclid is expected to 
have similar performances l lLaureiis et al. 2011). 

We stress that the resulting S/N ratio values should be deemed 
accurate only for multipoles up to ^ ~ 2000 — 3000. For angu- 
lar scales smaller than that non-Gaussian errors due to the non- 
linear growth of structures, which have not been taken into account 



4c* Jo \fKix)'^J a^ix) 



(32) 



(33) 



where the integration kernels are now defined as 

W.ix)= rn.ix')i^^^dx'. 

While previously the redshift distribution was normalized to unity 
over the complete redshift range, now we must normalize to unity 
the redshift distribution in each redshift bin, so that 



m{x)dx=l- (34) 
We did not attempt here a full tomographic analysis, since the 
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resulting gain in discriminatory power is likely not enough to dis- 
tinguish models 4 and 6, while the remaining models should be 
relatively easy to distinguish by simply using the weak lensing 
power spectrum alone. We did however compute the power spec- 
tra resulting from three different source redshift bins, in order to 
verify whether one of them would give a markedly stronger signal 
than the others in order to focus observational efforts on that red- 
shift range. The bins adopted are [0, 0.5], [0.5, 1], and [1, +oo]. In 
Figure [T5] we show the ratio of the power spectra restricted to a 
specific source redshift bin computed for each of the six oscillat- 
ing quintessence models in this work to the corresponding quantity 
evaluated in the framework of the ACDM cosmology. As can be 
seen, the qualitative behavior is the same as for the full effective 
convergence power spectrum, although some quantitative differ- 
ences exist. Specifically, the impact of oscillating dark energy is 
somewhat larger for the high redshift bins, having z > 0.5. How- 
ever it is a relatively small effect, changing the deviations with re- 
spect to the ACDM cosmology of a few percent at most. 



5 CONCLUSIONS 

In this work we considered structure growth in six different dark 
energy models characterized by an oscillating equation of state pa- 
rameter While many authors studied the expansion history of 
the Universe implied by these models, here we performed one step 
further by investigating the consequences of such models on the 
formation of non-linear structures. The main idea was to explore 
cosmological probes potentially capable of distinguishing one or 
more of these models from the concordance ACDM cosmology. 
To that aim, we studied several observables, ranging from the ISW 
effect to the cluster mass function and the cosmic shear power spec- 
trum. Our main conclusions can be summarized as follows. 

• No cosmologically relevant quantity shows oscillations as a 
function of redshift for any of the models considered in this work. 
This is a consequence of the fact that observables are given by in- 
tegrals over w[z), so that any feature in the latter function is effi- 
ciently smoothed out. There are very slight hints of a wiggle only 
in the ISW effect, in the age of the Universe, and in the deceleration 
parameter, at a level that is however likely impossible to detect. 

• We estimated the redshift drift, that is the variation of the cos- 
mological redshift of a source due to the expansion of the Universe 
in the various models considered here. The impact of oscillating 
dark energy can reach up to ~ 10% at 2 ~ 1 for the most extreme 
cosmologies, and stays at the level of a few percent all the way up 
io z ^ 10. At these high redshifts the impact of peculiar motions 



is highly negligible, thus allowing to probe the expansion history 
with very high accuracy. 

• The critical linear density contrast for spherical collapse 5c(z) 
and the virial overdensity Av(2) are always quite similar to the 
corresponding quantities evaluated within the fiducial cosmology, 
with deviations being at the level of a few percent for z < 2 — 3. 

• As naively expected, the growth factor of models having a 
larger amount of dark energy at early times is slightly larger than 
for the ACDM cosmology, and vice versa. Differences are at the 
level of a few percent also in this case. 

• The impact of oscillating quintessence on the mass function 
increases with both mass and redshift, reaching ~ 30% or more 
for M > lO^^Mo/ft at 2 = 2. Interestingly, models with a high 
frequency of oscillations show an increment in the abundance of 
cosmic structures, while models with a more regular evolution of 
w[z) show a decrement. This is because in models with high fre- 
quency of oscillations the amount of dark energy is higher than for 
the cosmological constant case and structures need to growth faster 
to compensate it (see Figures |2l|6] and [8}. The resulting effect on 
the redshift distribution of cluster catalogs depend heavily on the 
selection function, however it is at a level likely to be detectable 
with future wide cluster surveys. 

• Cosmic shear is affected at the level of ~ 10 — 15% at inter- 
mediate/small angular scales. Given the sky coverage, sensitivity, 
and PSF stability of future wide field weak lensing surveys, at least 
some oscillating dark energy models will be discriminated form 
the concordance ACDM cosmology by using the power spectrum 
of effective convergence alone. 

It is also worth to study what happens to our analysis if the 
dark energy equation of state is still oscillating, but is decreasing 
when the redshift increases in the vicinity of z = 0, opposite to 
what we assumed now. In order to achieve that we added a phase 
6 = n. We limited our new analysis only to the first three models 
and we consider the evolution in time of the growth factor, of the 
overdensities Sc and Av and the mass function at the four redshift 
considered, namely z = 0, 0.5, 1, 2. We observed that for all the 
quantities considered, the behaviour is qualitatively the same but 
the effect is much smaller than before. We can therefore conclude 
that for the oscillating dark energy models, a first oscillation in the 
equation of state with decreasing values of the equation-of-state pa- 
rameter makes the effects of dark energy smaller than before. This 
is due to the fact that the amount of dark energy in this situation is 
smaller than before, therefore structures do not need to grow as fast 
as it was before in order to compensate for it. 

It is interesting t o com pare our results with those of 
iMignone 8i BartelmannI ( |2008|) . In that work the authors showed 
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that even if the expansion rate of the Universe has a sudden tran- 
sition (that would require an even stronger transition in the dark 
energy equation of state parameter w{z)) this would hardly show 
up in cosmological observables such as the luminosity distance of 
SNe la. This agrees with our findings, namely since actual observ- 
ables are given by at least a double integral over the function w{z), 
any feature of the latter is easily smoothed out. A partial exception 
to this is the deceleration parameter. This is the second derivative 
of the expansion factor, hence it does retain some of the oscillatory 
behavior of the dark energy equation of state parameter 

We conclude by noting that, although it is virtually not pos- 
sible to find traces of oscillatory behavior in cosmological observ- 
ables, oscillating dark energy models do induce some specific mod- 
ifications in the number counts of massive clusters and the power 
spectrum of cosmic shear that will likely be detectable by future 
cosmological surveys. This paper hence gives an additional con- 
tribution to the study of observational signatures of a d3mamically 
evolving dark energy component, which is a fundamental field of 
study in order to better understand the nature of quintessence itself. 
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APPENDIX A: CODE DETAILS AND TESTING 

In this Appendix we explain in detail how we solved the Eqs. ((Sjl 
and (|6} for the non-linear and linear evolution of matter density 
fluctuations. Because all the terms involving the dark energy equa- 
tion of state present oscillating behavior, particular care has to be 
taken in order to achieve numerical convergence in the results. 
Since Eq. (|6]l is an ordinary second order linear differential equa- 
tion, in order to solve it we need to provide the initial conditions 5i 
and 5[. We assume that at early times the solution is a power law 
Si = a" (implying 5[ — nSi/a) and we insert this ansatz into the 
differential equation. By evaluating it at the equivalence scale fac- 
tor acq, we obtain a second order algebraic equation for n which 
can be easily solved. Since at early times all the models are very 
well approximated by an EdS model, n differs from unity only by a 
few percent, more so for models having more dark energy at early 
times. We parametrized the initial density contrast as a power law 
only in order to easily determine the velocity. 

To solve Eq. (|5} we have to provide as before two initial condi- 
tions, one for 5 and one for S' . At early times, we can safely assume 
that the evolution of density perturbations is linear and the solution 
can be written in the same form as done before. As explained in the 
main text, we need to determine the initial overdensity making the 
perturbation diverge at a given time. To do so we fix the time (scale 
factor) corresponding to the collapse and run a root-search algo- 
rithm. At each step Si and 5{ are automatically updated. The values 
corresponding to the divergent result of the non-linear evolution 
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are then used as initial conditions for the linear equation in order 
to determine the l inear overdensity parameter 5c- Despite the fact 
that, compared to'P ace et alJ (|20lrf), we do not have a fixed value 
for the initial velocity of the perturbations, we could reproduce the 
results presented there, implying that the role of the initial velocity 
is rather marginal. The advantages of doing so are twofold: on one 
side the new velocities are formally correct and can vary according 
to the cosmological model and on the other side the code is numer- 
ically much more stable. One instability problem is due to the fact 
that formally a divergent value needs to be infinite and this can not 
obviously be satisfied from the numerical point of view. We there- 
fore assume that divergence occurs when S 10^. This makes Sc 
artificially increase with the redshift also for an EdS Universe (see 
Figure[8ll. We checked that having the initial velocity related to the 
initial overdensity makes this problem much less severe. We run 
the code to determine the value of 5c up to Zc — 50 and we saw 
that first of all the increase is very mild, and second the numerical 
value is higher than the analytical one by only 0.5% at Zc — 50. 
For Eq. (O the initial scale factor is ai = 5 x 10~^. 

A crucial point for having numerical convergence of the re- 
sults is to perform an accurate numerical integration for the equa- 
tion determining the time evolution of the dark energy density. 
To perform the integral in Eq.[T2l we made sure to have a numerical 
accuracy of at least 10~^ when compared with the analytical result. 
We also checked that different integration methods would give the 
same result as indeed was the case. 



